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Abstract 

We propose a new method to measure the system noise temperature, T sys , using 
a 2-bit analog-to-digital converter (ADC). The statistics of the digitized signal in a 
four-level quantization brings us information about the bias voltage and the variance, 
which reflects the power of the input signal. Comparison of the variances in hot and 
sky circumstances yields T sys without a power meter. We performed test experiments 
using the Kagoshima 6-m radio telescope and a 2-bit ADC to verify this method. 
Linearity in the power-variance relation was better than 99% within the dynamic 
range of 10 dB. Digitally measured T sys coincided with that of conventional measure- 
ment with a power meter in 1.8-% difference or less for elevations of 10° — 88°. No 
significant impact was found by the bias voltages within the range between —3.7 and 
+ 12.8% with respect to the threshold voltage. The proposed method is available for 
existing interferometers that have a multi-level ADC, and release us from troubles 
caused by power meters. 
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1. Introduction 

System noise temperature (T sys ) is an important parameter for amplitude calibration in 
radio astronomy observations. This value represents the total equivalent system noise power, 
measured in the unit of Kelvin, including atmospheric radiation and attenuation, antenna noise, 
spill over, amplifier noise, loss in transmission, etc. For most radio sources and receivers, the 
received signal is dominated by the system noise power. Calibration for T sys is necessary to 
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measure the flux density or the antenna temperature of radio sources from the received signal. 

A conventional T sys measurement has been carried out using an analog power meter. One 
can estimate T sys by comparing the output signal powers between two or more circumstances 
with different input signals whose power is known. When the telescope points toward off-source 
blank sky, the output power, W^ky, is given as W s k y = kv,GBT SJS) where fee is the Boltzmann 
constant, G is the gain of the receiving system, and B is the bandwidth. In the case that the 
receiving system is covered with an absorber at the temperature Thot, i-e. "hot load", the output 
power, Whot, will be Whot = k-QGB{T SJS + Thot). Here Thot is the absorber temperature that is 
approximately equal to the room temperature of ~ 293 K. Using these two measurements we 
can estimate T sys as 

rn _ ^sky rp / 1 \ 

sys " w hot - w sky Ihot - [i) 

Precision, stability, and dynamic range are crucially required for the power meter to 
attain accurate T sys . The accuracy in T sys depends on the precision of the power meter. Since 
Whot is usually several times greater than PU s ky, the power meter must keep linearity in a power 
range as wide as several dB. A small non-linear characteristic of the power meter can result in 
systematic error in T sys . 

It is troublesome to maintain high precision and linearity in an analog sensor eternally, 
especially for radio telescopes in inaccessible environments such as space orbits, polar regions, 
and high-altitude plateaus. If a power meter device can be eliminated by an alternative method 
to measure T sys , we can reduce costs, troubles, maintenance frequency, etc. 

We propose a simple method to measure T sys using an analog-to-digital converter (ADC), 
instead of a power meter. Most of modern interferometers and radio telescopes are equipped 
with a high-speed ADC in each antenna element to acquire received signal in digital form that is 
processed in correlators and spectrometers. Utilization of digitized signal requires no additional 
costs. Since the digitized signal responds to voltage of the analog signal, its statistics represents 
the mean and variance (i.e. power) of the input signal. If the statistics reflects sufficiently 
accurate and reliable power in a sufficiently wide power range, ADCs can be the alternative. 

ADCs for radio interferometers sample the input analog voltage and quantize into a 
few — two, three, four, or more — levels. Distribution of samples in these levels holds statistical 
information such as mean and variance of the input signal. Coarse quantization can affect the 
performance in power measurements; two-level quantization is obviously degenerate in ampli- 
tude of the analog voltage and is useless in T sys measurements. Four-level (2-bit) quantization 
holds some information about amplitude. This quantization is commonly applied in radio in- 
terferometers such as ALMA (Okumura et al. 2002) and CARMA (Rauch et al. 2006), and 
VLBI such as VERA (Kobayashi et al. 2003), JVN (Doi et al. 2006), VSOP (Hirabayashi et al. 
1998), and VLB A (Romney 1988). 

In this paper we propose a method of measuring T sys using a four-level (2-bit) ADC. 
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Realistic performance of the method, tested with the Kagoshima 6-m radio telescope (Omodaka 
et al. 1994), is also reported. 



2. Method 

Thresholds in four-level quantization of ADCs are set at the voltages of ±Vo and 0, and 
the input analog voltages are encoded into 00 , 01 , 10 , and 11 as shown in figure 1 (Thompson 
et al. 2001). Distribution of samples in these codes reflects power and bias of input signals. 

We assume that the input signal is white noise whose probability density function of 
the voltage, p{V), follows the normal distribution with the mean (bias) voltage of \i and the 
standard deviation of a as 

~2a~ 2 

This assumption is justified for white noise that has a flat spectrum across whole bandwidth 
(Papoulis and Pillai 2002). The probability in each digitized code will be 
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These probabilities are constrained as P o + Poi + P 10 + Pa — 1- The standard deviation, a, and 
the bias voltage, /i, can be estimated by using the three independent observable values of Poo, 
P t, and P i0 . These probabilities are estimated by statistics of quantized signal with an ADC 
as P m i = n ml /n, when large enough number of samples are accumulated in each quantization 
level. Here, n ml is the number of samples in the quantization level of ml (m and 1 stand for the 
most and least significant bits, respectively) and n is the total number of samples. Non-linear 
least squares analysis works to obtain the parameters of and ai from the ADC statistics of 
P 00 , Poi and P 10 . 

The power of the input signal is proportional to the variance, a 2 , as a 2 = aW, where a 
is a proportional constant. Thus, equation 1 will be expanded as 
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Here, we denote Thot as a system noise temperature that is measured by the ADC statistics. 
Equation 4 indicated that Thot is calculated by ao, that can be derived from ADC statistics, in 
hot and sky circumstances. 

3. Tests and Results 

We examined the new T sys measurement by comparing with the conventional method 
using a power meter in four aspects — linearity, precision, dynamic range and robustness against 
input bias voltage. Linearity is estimated by relationship between the variance, a 2 derived from 
the ADC statistics and the power of input signal measured by a power meter. The proportional 
coefficient, a, should be constant within a wide range of input power in an ideal linear system. 
Dynamic range is evaluated by the width of the range where the system keeps proper linearity. 
Precision of digital T sys is evaluated by departures from the conventional T sys value measured 
simultaneously in the same configuration. Robustness against input voltage bias is tested by 
variation of a 2 obtained from bit distribution of ADC, for various input bias voltage in ADC. 

We used the Kagoshima 6-m radio telescope for these test experiments. The configura- 
tion is shown in figure 2. The antenna pointed to blank sky at various elevations to generate 
system noise. The feed horn inside the antenna was covered with a hot load when we measured 
PVhot- Received signal at 22.2 GHz was downconverted to intermediate frequency (IF) with the 
1st and 2nd local oscillators at 16.8 and 5 GHz, respectively. A step attenuator was used to 
adjust the IF power at a desired level. The IF signal at 256 - 512 MHz was bifurcated into 
two streams. One signal was digitized by the analog-to-digital converter (ADC), ADS-1000 
(Nakajima et al. 2001), at a sampling frequency of 512 Msps with four-level quantization. The 
other was used to measure the power in analog form with the power sensor, HP 8481 A attached 
to the power meter, HP 437B. 

Four experiments were carried out to test the aspects mentioned above. 

3.1. Distribution of codes 

We examined statistics of samples coded in 00, 01, 10, and 11 in various input power 
and bias offset. The total number of samples was 512 x 10 6 in 1-sec accumulation. Figure 3 
shows distributions of acquired samples, as a function of input power, without bias and with 
the maximum bias offset of h/Vq = 0.13. 

The distribution was symmetric without bias and became significantly asymmetric with 
the bias. Percentages of upper levels of 00 and 11 were 0.00949 % (48570 samples) and 0.00975 % 
(49925 samples) without a bias and 0.00097 % (4963 samples) and 0.04864 % (249027 samples) 
with the maximum bias offset. These shares increased at higher input power, as expected in 
equations 3. 
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3.2. Power and digital variance 

To verify the linearity between the input analog power and the digital variance, a 2 , we 
compared them in various input power range by using a step attenuator between the bandpass 
filter (BPF) and the power splitter to adjust input signal power. The input signal power was 
varied in the range from 0.03 mW (—15 dBm) through 7.9 mW (+9 dBm) with a 1-dB interval 
for the zero bias offset case. 

For the case with a bias offset of /z/Vo = 0.13, the maximum input power was capped at 
2.8 mW (4.7 dB) to avoid saturation. We realized that this limitation was far from saturation 
and was not necessary after the measurement. However, the linearity in the biased case was 
grasped with the capped power range. 

Figure 4a shows the relation between the analog power and the digital variance normal- 
ized by the threshold as (a/V ) 2 . The best-fit linear regressions of (cr/Vo) 2 = a/V 2 W yielded 
a/V 2 = 2.548±0.001 and a/V 2 = 2.524±0.0007 for bias offsets of fi/Vo = and 0.13. Departure 
from the linear regression was 4.0% and 4.6% at the minimum and maximum input power of 
0.03 and 7.9 mW in the case of zero bias offset. The departure in the biased case is not signifi- 
cantly different from the unbiased case. The departures remain < 1.0% in terms of peak-to-peak 
in the input power range of 0.09 and 1.8 mW, that yields a dynamic range of 13 dB. 

3.3. Bias endurance 

An additional linearity test was taken under a bias voltage in the input analog signal. 
Instead of adding a bias voltage in the input analog signal, we adjusted the threshold voltages 
of the ADC. The relative bias shift, fi/Vo, were estimated by the statistics of digitized signals. 

Figure 46 shows the linearity relation, as same as figure 4a, with the bias shift of h/Vq = 
0.13. In this test the input power was changed from —40 through —20 dBm. The best-fit linear 
regression was aW/V 2 = (1 ± 0.0007) {a/V ) 2 . 

The ratios, a, in various bias voltages are plotted in figure 5. No significant departure 
greater than 0.4% from 1 was found in the tested relative bias range of —0.04 < fi/Vo < 0.13. 

3.4- Analog and digital T sys measurements 

We examined accuracy of digitally-estimated T sys by comparing with T sys of the conven- 
tional method, in various elevations. Whot was measured once at the elevation of 88°, that was 
the upper limit of the 6-m telescope. W s ^y was taken at 88° and from 85° through 10° by a 5° 
step. Statistics of digital samples was also acquired at each elevation. 

Figure 6 shows the comparison between T sys and T sys . The linear regression resulted in 
^sys = (0.988 ± 0.019)T sys + 0.3 ± 6.6. Departure of T sys /T sys from unity was less than 1.9% in 
whole T sys range from 303.7 through 1150.4 K. 
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4. Discussion 



4-1. Performance in system noise measurements 

As shown in figure 4 a 2 proportionally responses to the input power. Departure from 
the ideal linear response is 4% at the lowest input power and 4.6% at the maximum power of 
7.85 mW (+9 dBm). The departure from the diagonal line is within 1% (between —0.55 and 
— 1.55 %) if we use the input power within 0.09 mW to 1.8 mW. Thus, the linearity is as good 
as 99% in this dynamic range of 13 dB without any correction. This error range is comparable 
with the accuracy of the power meter, specified as ±0.5%. 

The departure increases as input power exceeds 1.8 mW. This behavior can be caused 
by saturation in the power meter. Relevant calibration curve would reduce the departure and 
overcome the nonlinearity. 

Since the power ratio between hot and sky inputs is usually less than 10 dB, the linearity 
without correction is enough for Tsys measurements with the accuracy of 1%. 

The dynamic range of 13 dB for the 99% linearity is acceptable for usual Tsys measure- 
ments using hot and cold loads. For lower-noise system, wider dynamic range is needed. When 
Tsys is better than 32 K, hot-sky ratio exceeds 10 dB. Relevant calibration curve would expand 
the dynamic range and allows us to measure such low-noise system. 

4-2. Accuracy in system noise measurements 

Error in T sys is estimated from equation 1 as 



ST sys \ 2 _ / W sky \ 2 



T sys J \W hot -W skYi 



'SW aky \ 2 , [5W hot 



w sky J ' \w hot 1 1 (5) 

where 5W sky and 5Whot are errors of power in sky and hot circumstances. In the case of digital 
T sys measurements, the error in power consists of error in a 2 and variance in a within certain 
range of hot and sky powers. 

Statistical ambiguity in bit distribution causes the error in a 2 . Uncertainty in n ml is 
estimated by variance of the binominal distribution, nP ml (l — P m i), so that the error in P m ± is 
given by 

3 -Pml(l - -Pml) /„n 

= • (6) 

The number of sample, n, equals to the product of the sampling rate and the integration time. 
In our experiment the sampling rate of 512 Msps and the integration of 1 second brought 
n = 5.12 x 10 8 and resulted in a% ~ 10~ 9 . 

■"nil 

Variance in a is estimated by the results of the linearity test as shown in figure 4, that 
is less than 1% within the range of 10 dB. This term is much larger than the error in a 2 , hence, 
is the major determinant in 5W sky and <5PUhot unless P ml < 10 -5 . Using equation 5 the total 
error in digital T sys measurement is estimated to be ~ 1%. As shown in figure 6, T sys measure 
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by the bit distribution coincides the conventional measurements within 1.88%. The maximum 
error appears at the lowest elevation of 10° where the factor W s k y / (Whot — Wsky) is largest. The 
magnitude of departure is reasonable as estimated in equation 5, with some additional error 
caused by the power meter and different signal path after the power splitter. 

4-3. Robustness against bias 

The power ratio differs only 0.2% when added bias of —0.037 to +0.128 Vq. The results 
simply show that new digital power measurements is robust enough against input voltage bias. 
As shown in equation 3 a and \i can be well orthogonally determined by the non-linear least 
squares. This is why T sys was not significantly affected by the bias. Larger bias, ^> \Vq\, 
however, can cause very small probability in P 00 or P n that can result in relatively large error 
in T sys . 

4-4- Inconvenience as a power meter 

Although the new method works for the T sys measurement, it is not suitable as a general- 
purpose power meter. The input signal voltage is presumed to follow the normal distribution 
in the estimation of a. The assumption is improper unless the input signal is a white noise, 
such as the system noise. A monochromatic non-biased signal at a particular frequency with 
the amplitude of Vo, for instance, gives a probability density function as p(V) = ^{V 2 — V 2 ) 1 / 2 , 
that is different from the normal distribution. A 4-level ADC with the threshold voltage of Vq 
will exhibits a non-linear response to such a signal and the bit distribution obviously freezes 
when Vq < Vq. The bit distribution is effective only when the probability distribution is known 
and spans the range of digital states. 

5. Summary 

We examined performance of T sys measurement using a 2-bit ADC, and verified that the 
new method is available as an alternative of conventional measurement using a power meter. 
The statistics of digitized samples showed a sufficiently linear response whose departure from 
power-meter measurement was less than 1% in a dynamic range of 13 dB. Although we added 
bias offsets in input voltage up to 0.128 Vq, the linear response was kept in the same dynamic 
range. The linearity and dynamic range yielded accurate T sys measurements that differed < 1.9% 
from conventional measurements using a power meter at various antenna elevation angles. 

Modern radio telescopes and interferometers that are equipped with a high-speed ADC 
with 2-bit quantization can use this new method. ADCs with 3- or more-bit quantization are 
expected to perform better accuracy and a wider dynamic range with an appropriate calcu- 
lation of bias offset and standard deviation similar to this method. This method is available 
without any additional hardware. Since this method counts digitized signals that will be used 
in spectroscopy or interferometry, it would be less affected by systematic errors that is cased 
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by non-linearity and different passband characteristics in a different signal path to a power 
meter used in conventional T sys measurements. Thus, this method can be a new standard of 
T sys measurements in radio astronomy. 

We thank Dr. J. D. Romney who reviewed this paper and gave positive suggestions to im- 
prove the manuscript. The Kagoshima 6-m telescope is supported by the National Astronomical 
Observatory of Japan and operated by staffs and students of the Kagoshima University. 
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Fig. 1. Four-level encoding and probability density distribution. Horizontal axis indicates the input 
(analog) voltage. Top diagram shows the relation between the input voltage and the output codes 00 , 
01, 10, and 11, with the thresholds of — Vo, 0, and Vo- Bottom diagram shows the probability density 
function, assuming the normal distribution with the mean and standard deviation of \x and a. Probabilities 
of the output codes corresponds to the filled area in each level. 




Fig. 2. Test configuration. The Kagoshima 6-m telescope pointed various elevations of the blank sky. 
Received signal was amplified in the low-noise amplifier (LNA) and was downconverted with the 1st and 
2nd local oscillators (LOs) at 16.8 and 5 GHz to an IF signal. The IF signal was filtered in the bandpass 
filter (BPF) and bifurcated by the power splitter. It was quantized by the analog-to-digital converter 
(ADC). The IF power was measured by the power meter (PM) in analog form at the same time. The 
hot load can cover a feed horn inside the telescope. The IF power was adjusted at a desired level by the 
attenuator (ATT). 
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Fig. 3. (a) Distribution of acquired samples coded as 00, 01, 10, and 11 without bias offset. 
Shares in the outer codes of 00 and 11 were 0.00949 % and 0.00975 %, respectively, at the 
lowest input power of 0.07 mW and increased at higher input power. (b) Distribution with a 
bias of n/Vo = 0.13. The distribution became significantly asymmetric across the codes. Shares 
in 00 and tt 11 were 0.00097 % (i.e. 4963 samples in 1 sec) and 0.04864 %, respectively. 
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Linearity of Power Measurements Linearity of Power Measurements 
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Fig. 4. (a) Linearity without offset bias. (Upper Panel): Relation between the analog input power, 
W [mW], and the variance normalized by the threshold voltage, (a/Vo) 2 , calculated from the statis- 
tics of digitized signals. The solid line indicates the best-fit linear regression of (er/Vb) 2 = cvW/Vq, 
where a/Vg = 2.548 ±0.0001 mW" 1 . (Bottom): Ratio of the variance to the linear regression. The 
vertical axis indicates o~ 2 /(aW). The departure from the linear regression becomes larger, up to 
4.0% and 4.6% at the lowest and the highest input power, respectively. It remains almost con- 
stant with the peak-to-peak range of 1.0% in the input power range between 0.09 and 1.8 mW, 
that yields a dynamic range of 13 dB. (b) Linearity with a bias offset of fi/Vo = 0.13. The lin- 
ear regression shows a = 2.524 ± 0.0007. The departure from the linear regression remains < 1.0 % 
(p-p) in the input power range between 0.09 and 2.2 mW that yields a dynamic range of 14 dB. 
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Fig. 5. Bias endurance of power measurements. The variance-to-power ratio 

(a/Vo) 2 /aW was measured at several bias offset between —0.04 < (j,/Vq < 0.128. 

Throughout bias offset, the ratio was almost unity with no errors beyond 0.4%. 
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Fig. 6. (Top:) Comparison of digitally measured T sys with the conventional T sys using a power meter. 
T sys ranged between 303.7 K and 1150.4 K for elevation angles of 88° and 5° of the Kagoshima 6-m radio 
telescope. (Bottom:) The ratio of digital T sys to conventional T sys measurements. Relative error from 
unity was less than 1.9%. 
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